m  ri7 


RIA-80-U302 


USADACS  T»chntc*l  Ubrtry 


5  0712  01013444  2 


TECHNICAL 
t  LIBRARY 


COPY  NO. 


TECHNICAL  MEMORANDUM  1891 

RELY, 

A  PROGRAM  FOR  ESTIMATING 
OVERALL  SYSTEM  RELIABILITY  BASED  ON 
COMPONENT,  SYSTEM  AND 
FLIGHT  DATA 

THOMAS  A.  NEEF 

MARCH  1969 


THIS  DOCUMENT  HAS  BEEN  APPROVED  FOR  PUBLIC 
RELEASE  AND  SALE;  ITS  DISTRIBUTION  IS  UNLIMITED 


PICATINNY  ARSENAL 

DOVER,  NEW  JERSEY 


TECHNICAL  MEMORANDUM  1891 


/ 


RELY, 

A  PROGRAM  FOR  ESTIMATING 
OVERALL  SYSTEM  RELIABILITY  BASED  ON 
COMPONENT,  SYSTEM  AND 
FLIGHT  DATA 


BY 

THOMAS  A.  NEEF 
MARCH  1 969 


DATA  PROCESSING  SYSTEMS  OFFICE 
PICATINNY  ARSENAL 
DOVER,  NEW  JERSEY 


This  document  has  been  approved  for  public  release  and  sale;  its 
distribution  is  unlimited 


TABLE  OF  CONTENTS 


PAGE 

ABSTRACT  2 
INTRODUCTION  3 
GENERAL  PROGRAM  OUTLINE  h 
COMPUTER  FROGRAM  DESCRIPTION  11 

/ 

(INCLUDING  DESCRIPTION  OF  SUBHIOGRAMS)  12 
PROGRAM  OPERATION  18 
SAMPLE  CASE  (INCLUDING  SAMPLE  INPUT/OUTPUT)  21 
APPENDIX  I  -  PROGRAM  LISTING 


T 


1 


I 


ABSTRACT 


This  report  documents  a  program  for  estimating  overall  system 
reliability  by  judiciously  combining  individual  component  data, 
laboratory  test  data  and  finally  flight  data  as  gathered  at  two 
different  periods  in  time.  All  data  used  in  arriving  at  this  single 
reliability  estimate  is  assumed  to  be  attribute  in  nature. 


j 
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INTRODUCTION 


In  July  1968,  the  Munitions  Command  requested  the  Mathematical 
Analysis  Division  of  the  Data  Processing  Systems  Office  to  develop  a 
computer  program  (using  their  prescribed  methodology)  to  estimate  over¬ 
all  system  reliability.  Such  a  program  has  been  created,  herein  called 
RELY,  which  requires  as  input  six  sets  of  distinct  attribute  type  data 
which  should  commonly  arise  during  the  development  phase  of  a  missile 
assembly  or  subassembly  program.  The  six  sets  of  basic  inputs  consist 
of  the  three  types  of  data  listed  below  gathered  at  two  distinct  time 
periods: 

(1)  Component  data  as  naturally  arising  from  individual  component 
tests. 

(2)  Laboratory  systems  data  arising  from  testing  a  component 
assembly  (i.e.  circuit)  in  a  laboratory  and 

(3)  Flight  data  consisting  of  data  resulting  from  actual  test 
firings  containing  all  circuitry. 

The  program  judiciously  combines  this  data  so  as  td  produce  a 
single  estimate  for  overall  system  reliability. 

This  report  describes  in  detail  this  computer  program  and  includes 
a  complete  description  of  input-output  formats,  a  discussion  of  the 
program  logic  and  a  sample  case.  It  should  be  emphasized  however,  that 
this  report  does  not  describe  the  foundations  of  the  methodology 
utilized.  For  this  the  interested  reader  is  referred  to  Reference  l.» 
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GENERAL  PROGRAM  OUTLINE 


To  simplify  the  discussion  the  following  mathematical  notation  has 
been  adopted. 

Mathematical  Notation: 

A.  .  A  =  N  refers  to  number  of  tests  conducted  at  i  th 

ij  — 

stage  at  time  T^ 

(basic  inputs)  A  =  C  refers  to  number  of  successes  resulting  from 
testing  at  the  i  th  stage  at  time  T^ 

Subscript  definition: 


i  =  c  (Component  stage)  refers  to  component  data  . 
i  =  s  (£ystem  stage)  refers  to  system  data,  i.e. 
resulting  from  testing  circuits  comprised  of  the 
above  components 

i  =  f  (flight  stage)  refers  to  flight  data,  i.e. 
resulting  from  flight  tests  containing  the  circuitry 
above 

j  =  o  refers  to  time  Tq  ' 

j=l  refers  to  time  T  where  T  <  T 

1  o  i 

Thus  N  refers  to  the  number  of  systems  tests  conducted  at  time  T  , 

Sl  1 

while  refers  to  the  number  of  flight  successes  obtained  at  time  Tq. 
K.  refers  to  a  weighting  factor  that  weights  the  less 

significant  Tq  data  when  combining  with  T  data,  i 
again  refers  to  component  (i=c),  systems  (i=s)  or 
flight  (i=f)  data. 


k 


A. 

leq 


K. 

leq 


refers  to  the  equivalent  number  of  tests  (A=N:, 

or  equivalent  number  of  successes  (A=C)  for  either 

component,  system  or  flight  data,  (i=c,  s,  f 

respectively).  These  values  arise  by  appropriate 

combination  of  N  ,  C  ,  N  ,  C  ,  K  to  obtain  N 

coJ  co*  Ci*  ci’  c  ceq 

and  Cce(^,  or  more  generally  by  combination  of  A^q, 

A^  and  K^,  to  produce  A^  . 

refers  to  the  weight  given  to  the  equivalent  compon¬ 
ent  data,  when  combining  with  the  equivalent  system 
data  (i=s),  to  produce  combined  component-system  data. 
When  i=f,  this  variable  refers  to  the  combination  of 
the  component-system  data  with  the  equivalent  flight 
data  to  produce  the  overall  equivalent  number  of 
tests  and  overall  equivalent  number  of  successes 

C  used  to  arrive  at  P,  .  . 
eq  best 


best 


refers  to  the  overall  system  reliability  =  C  /N 

eq  eq 


The  interrelationships  existing  between  the  various  data  sets  are 
exhibited  in  Figure  1.  Figure  2  is  a  schematic  diagram  showing  the 
relative  significance  of  the  involved  data  sets. 

The  underlying  approach  to  arrive  at  P^es^.>  ‘the  single  estimate' 
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of  overall  reliability,  is  to  combine  all  like  data  sets,  degrading 
however,  the  less  significant  Tq  data,  when  combining  with  the  more 
significant  T ^  data.  Similar  degradation  factors  apply  when  combining 
the  less  significant  equivalent  component  data,  with  the  more  significant 
equivalent  system  data  and  the  more  significant,  yet,  flight  data. 

Each  such  combination  of  data  will  be  described,  thus  yielding  in 
the  process  a  general  program  outline. 


I.  Combination  of  component  data  N  ,  C  ,  N  ,  C  to  produce  N  , 

c  coJ  co’  Ci*  Ci  *  ceq’ 

C  . 
ceq 

A.  One  component  case. 


Here  N  ,  C  ,  N  ,  C  embocty-  the  inputs  to  this  calculation.  The 

CO  CO  c  i  c  ^ 

values  for  the  above  are  obtained  from  actual  testing,  witnessing  for 

example  C  success  out  of  N  tests  at  time  T  . 

co  co  o 


FIGURE  1 
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It  is  assumed  that  P,  the  true  probability  of  component  function¬ 


ing  is  ^-distributed  with  Ncq,  Cco  as  parameters,  i.e. 


Po(p) 


r  (n  +  2) 

co 


r(c  + 1)  r(N  -  c  + 1) 

A  CO  A  CO  CO 


C  N  -C 


(i) 


with  a  similar  expression  (p)  for  components  tested  at  time  T^. 

The  Tq  data  is  degraded  try  the  area,  Kc,  common  to  the  two  ^  -  distrib¬ 
utions,  as  expressed  in  the  following  equations  that  yield  the  equival¬ 
ent  component  data. 


N 


ceq 


N 


co 


K  +  N  : 
c  Cl' 


C  =  C  .  K  +  C 
ceq  co  c  cj 


(2) 


is  interpreted  geometrically  in  the  Diagram  below 


£  I 


x 
i  \ 


4 

/ 


I,.. 


DIAGRAM  1 


Thus  if  N  =  N  and  C  =  C  the  two  distributions  would  coincide 

CO  Ci  CO  Ci 

and  the  degradation  factor  would  be  1  (i.e.  no  degradation). 

B.  Several  Component  case. 

When  several  different  component  types  are  involved  and  related 

try  a  specific  circuit  equation,  calculations  additional  to  those  above 

are  necessary.  In  this  case  it  is  required  to  compute  a  N  and  C 

J  mi-  ceq  ceq 

for  each  component  type  and  combine  the  results.  Let  NQ  (k)  and 
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CCeq(k)  represent  the  above  values  referred  to  the  k-th  component  type. 
An  equation  similar  to  equation  (l)  is  constructed  with  these  values 


i.e. 


P  (P,k) 


r  (Nceq(k)  ♦  2) 


r(cceq(k)  +  d  r(Mceq(k)  -  oceq(k)  +  d 


(3) 


where 


(  1  -  P  ) 


N _ (k)  -  C  (k) 

ceq 


C_(k) 

Q  =  P 

This  represents  the  probability  density  P  of  the  k-th  component  type 


in  the  light  of  both  T  and  T  data.  The  calculation  to  combine  each 

N  (k)  and  C  (k)  to  produce  N  and  C  then  proceeds  in  a  Monte 
ceq  ceq  coq  ceq 

Carlo  fashion  involving  the  circuit. 

Values  are  sampled  from  each  of  the  (V  distributions^  and  in¬ 
serted  into  the  specified  circuit  equation,  from  which  the  probability 

of  circuit  functioning,  R  for  these  values,  is  obtained.  It  is  import- 

s 

ant  to  note  that  if  a  given  component  type  appears  several  times  in  the 
circuit,  the  corresponding  A-  distribution  is  sampled  only  once  per 
simulation.  Performing  this  simulation  M  times  thus  yields  M  values 
for  the  probability  of  circuit  functioning.  The  assumption  is  made 

f 

that  these  M  values  are  ^'distributed.  Consequently  the  mean  u. 
and  the  variance  C~2  of  these  numbers  are  computed  and  set  equal  to 
the  mean  and  variance  of  the  ^-distribution. 

*This  is  accomplished  by  numerically  determining  the  cumulative 
distribution  and  finding  the  P  corresponding  to  RN  =  f  k,X)  dX  where 
RN  is  a  random  number  between  0  and  1.  The  integral  is  tabulated  as  a 
function  of  P,  so  that  for  given  RN,  linear  interpolation  applies  for 
ascertaining  P. 
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This  allows  for  the  two  unknown  parameters  in  this  distribution  to  be 


solved.  These  unknowns  correspond  in  fact  to  the  N  and  C  of 

ceq  esq 

Figure  1.  Mathematically,  this  procedure  is  equivalent  to  the  equations 


that  follow.  ^ 

ZRs(i) 

i  =  1 


M  2  *A  i 

^]Rs(i)  -  2  Aps(i) 


M 


<t2  -  LlL 


i  =  1 


M 


(h) 


N _ =  &  (1  -  /U, )  -3 

_2 


ceq 


°ceq  '  (l  *P-  > 


(5) 


II.  Combination  of  system  data  N  ,  C  ,  N  ,  C  to  produce  N  and 

J  so’  so’  si’  si  1  seq 

C  , 
seq’ 

This  combination  is  accomplished  analogously  to  equations  (l) 
where  the  degradation  factor  is  again  the  area  common  to  the  two 
distributions .  Thus 


i'(n  +  2) 

SO 


r(c  +  i)  r  (n  -  c  +i) 

SO  so  so 


C  N  -  C 

P  ^^(l  —  p) 


f-  (P) 


r(»s,+  2) 


1  r(c  + 1)  r  (n  -  c  +  i) 

Si  SI  Si 


C  N  -  C 

P  Sl(l  -  P)  31  Sl 


K-  Miu{  (a  (P),6  CP>>  dP 

S  ~  O  VO  ’  '  x 


(6) 


N  -  N  . 
seq  so 


K  + 
s  s 


C  -  C  . 
seq  so 


K  +  C 

s  Si 


< 
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III.  Combination  of  flight  data  N^,  C^,  C^  to  produce  , 


Cfeq* 


This  again  is  completely  analogous  to  II  with  Ngo, 


replaced  by  N^o,  . ..,  C^  respectively. 


IV.  Combination  of  equivalent  data  produced  in  I,  II,  and  III  to 

produce  N  ,  C 
^  /  eq’  eq 

The  combination  of  N  ,  C  with  N  ,  C  to  produce  N  , 

ceq’  ceq  seq’  seq  c,seq’ 

Q 

Jc , seq  and  the  combination  of  these  last  two  values  with  N „  ,  C  „ 

’  feq’  feq 

to  produce  N  ,  C  are  also  completely  analogous  to  II  with  the  A 

6Q  6Q  ceq 

replacing  the  A  and  the  A  replacing  the  A  in  the  first  case  and 
^  b  so  seq  Si 

the  A  ^  replacing  the  A  and  the  A„  replacing  the  A  in  bhe 
c  y  seq  so  xeq  si 

second  case.  Thus  the  last  equation  in  (6)  becomes  C  =  C 

n  c,seq  ceq 

K  +  C^  when  combining  components  and  system  data,  and  becomes 
seq  seq 

C  =  C  .  K„  +  C~  when  combining  component-systems  data  with 
eq  c  y  seq  xeq  xeq 

flight  data. 


V.  Calculation  of  overall  system  reliability,  P^est' 

C 

P  =  _ 54- 


best 


N 


eq 


(7) 


For  output  purposes,  the  following  two  variable  are  also  computed  by 
the  equations  that  follow 


VARIANCE  =  (C  +  l)  (N  - 
eq  eq 


C  +  1) 

_ 


Mu  =  C  +  1 

e£.. 


COMPUTER  PROGRAM  DESCRIPTION 


A.  Description  of  MAIN  program 

This  Section  describes  the  computer  program  that  evaluates  the 
mathematical  model.  It  is  believed  that  sufficient  information  is 
contained  herein  along  with  COMMENT  statements  in  the  program  listing 
to  enable  one  to  understand  and  possibly  modify  the  program.  Figure  3 
is  an  overall  flow  chart  giving  the  sequence  of  the  computations.  Not 
shown  in  the  figure  is  the  fact  that  much  of  the  routine  calculations 
such  as  the  tabulation  of  the  distribution,  numerical  integration 
etc.  are  relegated  to  subroutines.  These  subroutines  though  briefly 
referred  to  in  the  text  are  more  fully  described  at  the  end  of  this 
Section. 

All  inputs  are  read  in  the  MAIN  program,  thereupon  one  set  of 
component  data  is  selected,  being  the  values  for  N  and  C  at  time  Tq  and 
values  for  N  and  C  at  T  .  Using  these  four  values  (labelled  in  the 
program  as  NCO,  CCO,  Nd,  CC1  respectively)  as  input  arguements  the 
value  KA  (=Kc  of  Diagram  l)  is  ascertained  from  a  subprogram  K.  The 
Tq  data  is  then  degraded  by  this  factor  and  the  result  added  to  the  T 
data  to  yield  two  new  intermediate  values,  SUB1,  corresponding  to  the 
number  of  trials  and  SUB2  corresponding  to  the  number  of  successes  as  in 
equation  ( 2 )  (page  7).  Using  these  intermediate  values  as  parameters  a 

distribution  (corresponding  to  equation  3)  is  constructed  and 
tabulated  in  subroutine  BETAF.  The  tabulation  is  effected  at  equally' 
spaced  abscissas  depending  upon  the  value  of  DELTAP  which  should  lie 
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Figure  3 


in  the  interval  [.01,  .5]  .  It  is  from  this  density  function  that  the 
probability  of  component  functioning  is  sampled  and  subsequently  in¬ 
serted  into  the  circuit  equation.  To  accomplish  this,  the  cumulative 
of  this  distribution  is  tabulated  and  placed  in  an  array  A(l)  (sub¬ 
routine  DQSF),  a  random  number,  (RDM(l)),  selected  and  the  correspond¬ 
ing  P  value  determined  by  linear  interpolation  in  subprogram  INTERP. 

This  calculation  is  repeated  for  each  component  type,  and  all  results 
are  ultimately  substituted  into  the  circuit  equation  from  which  a 
probability  of  circuit  functioning,  (RS),  is  computed.  At  this  point 
a  running  count  of  these  RS  values  for  output  purposes  in  histogram 

form  is  maintained,  along  with  running  summations  of  the  RS  values  and 

2 

RS  values  for  use  in  calculating  the  mean  jx,  and  variance  j  *  of 

/ 

these  generated  numbers  (equations  i;).  At  the  end  of  the  simulation 

loop,  equations  (5)  are  used  to  produce  N  and  C  which  are  labelled 

ceq  C6Q 

as  NCEQ  and  CCEQ  in  the  program.  Remaining  calculations  involving 
system  and  flight  data,  and  all  combinations  thereof  are  analogous  to 
what  has  been  described  and  again  typically  involve  use  of  the  BETAF 
and  K  routines.  The  notation,  too,  in  remaining  segments  of  the  MAIN 
program  is  quite  suggestive  so  that  the  reader  will  find  for  example 
Nso,  Cggqj  and  *best  3X6  represented  in  the  program  as  NSO,  CSEQ  and 
PBEST. 

B.  Description  of  Subprograms  (in  alphabetical  order) 

1.  FUNCTION  AREAF  (B,j) 
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Common  DELTAP 


Description  of  arguments: 

B  -  one-dimensional  equidistantly  tabulated  function  the 
area  under  which  is  to  be  computed. 

J  -  the  number  of  values  in  B. 

DELTAP  -  interval  width 
I  Operation: 

Calculates  the  area  under  B  by  Simpson's  rule,  (there¬ 
fore  J  must  be  odd).  The  output  is  a  scalar  representing 
the  area  under  the  curve. 

2.  SUBROUTINE  BETAF  (N,  C,  BETA,  J,  PB) 

Common  DELTAP 

Description  of  arguments: 

Input  Arguments : 

N,C  -  input  parameters  for  the  beta- distribution  equation 
Output  Arguments: 

BETA  -  one  dimensional  beta- distribution  array  produced 

f 

in  the  subroutine 

J  -  number  of  values  in  array  (dependent  upon  DELTAP) 

PB  -  one  dimensional  array  of  P-values  used  in  generating 
the  beta  distribution 
DELTAP  -  interval  width 
Operation: 

f 

Produces  an  equidistantly  tabulated  beta-distribution 
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using  the  following  equation: 


(3  =  G  PC  (l  -  P)  N  "  C 


(9) 


where  G  is  a  constant, 


r (N  +  2) 


G  = 


(10) 


r(c  +  i)  r(N  -  c  + 1) 


and  is  computed  in  subroutine  GAMA,  P  varies  from  0.0  to 


1.0  in  increments  of  DELTAP 

In  order  to  avoid  problems  arising  from  zero  bases  and/or 
zero  exponents  in  the  exponentiation  routine,  the  equation 
is  evaluated  for  one  of  four  cases  as : 

1.  N  and  C  are  distinct  and  nonzero. 

2.  N  and  C  are  identical  and  nonzero. 

3.  N  is  nonzero  but  C  is  zero, 
li.  Both  N  and  C  are  zero. 


3.  SUBROUTINE  DQSF  (H,  Y,  Z,  NDIM) 

This  subprogram  is  part  of  the  IBM  S|ystem/360  SCIENTIFIC 
SUBROUTINE  PACKAGE  and  a  detailed  description  'of  its  argu¬ 
ments  and  operation  can  be  found  therein.  Briefly,  however, 
this  subroutine  performs  the  integration  of  an  equidistantly 
tabulated  function  by  Simpson’s  Pule.  It  computes  a  vegtor 
of  integral  values  Z^  for  a  table  of  functional  values 
Y^,  i=  1,  2,  ...,  N  given  at  equidistant  points  X^  =  a  +(i-l)h 


(i=  1,  2,  3,  ...,  N) 


Ik 


h.  FUNCTION  GAMA  (ARG1,  ARG2,  ARG3) 
Description  of  arguments: 


ARG1  -  corresponds  to  N+2  in  equation  (10)  of  subroutine 
BETAF 

ARG2  -  corresponds  to  C+l  in  equation  (10)  of  subroutine 
BETAF 

ARG3  -  corresponds  to  N-C+l  in  equation  (10)  of  subroutine 


BETAF 


Operation: 

The  IBM  -  supplied  subroutine  GAMMA  evaluates  the  GAMMA 
function  for  a  given  argument  provided  the  argument  is 
less  than  or  equal  to  57.  The  subroutine  GAMA  is  designed 
to  reduce  each  of  the  arguments  ARG1,  ARG2,  ARG3,  to  57, 
if  necessary,  before  calling  GAMMA  for  each  of  them. 
Furthermore  it  computes  G 


r(ARGl) 


(11) 


r (ARG2)  r( ARG3) 


for  use  in  subroutine  BETAF.  Output  is  scala!r  value  for  G. 


5.  SUBROUTINE  GAMMA  (XX,  GX,  IER) 

This  subroutine  is  part  of  the  SCIENTIFIC  SUBROUTINE  PACKAGE 
except  that  was  changed  to  double  precision  for  this  program. 
Briefly  the  operation  of  this  subroutine  is  to  evaluate  the 
gamma  function  for  a  given  value  of  XX,  where  V  (x)  is  defined 
for  X  >  0  by  • 
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which  satisfies  the  recurrence  relation  r'(X)  =  (X-l)r(X-l) 
which  defines  X*(x)  for  ary  X  non-negative  integer  i.e. 

r(x)  =  (x-D! 

6.  FUNCTION  INTERP  (i,  J,  F,  YFL) 

Common  /BLKL/AX 

Description  of  arguments: 

'  I  -  component  type  designation 

J  -  number  of  values  in  ITH  row  of  two-dimensional  array 
AX  corresponding  to  array  AREA  of  the  MAIN  program,  in 
which  the  scalar  input  YFL  will  be  interpolated. 

F  -  the  one- dimensional  array  from  whose  elements  the 
interpolated  value  is  to  be  calculated. 

YFL  -  number  whose  corresponding  value  is  to  be  found  by 
interpolation 

AX  -  subroutine  name  for  array  AREA  of  MAIN  program  which 
stores  cummulative  for  each  component  type. 

Operation: 

Straight-forward  linear  interpolation  is  effected  herein. 

7.  FUNCTION  K  (NO,  CO,  Nl,  Cl) 

Common  DELTAP 

Description  of  arguments 

NO,  CO  -  number  of  trials  and  number  of  successes  at  time 
Tq  (or  analogous  parameters) 
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Nl,  Cl  -  number  of  trials  and  number  of  successes  at  most 
recent  testing  T^  (or  analogous  parameters) 

DELTA P  -  interval  width 
Operation: 

K  is  designed  to  produce  a  proportionality  facotr  for  de¬ 
grading  the  less  significant  data,  NO,  CO  (before  it  is 
combined  with  the  more  significant  data  Nl,  Cl  in  MAIN 
program).  This  is  accomplished  by  first  tabulating  the 
beta-distributions  for  NO,  CO  and  for  Nl,  Cl  by  calling 
subroutine  BETAF.  The  area  under  the  intersection  of  these 
two  curves,  CMAKEA,  is  found  next  by  comparing  the  two  beta 
curves  point-for-point,  in  each  case  taking  the  smaller  of 
the  two  values  (ordinates),  putting  it  into  a  new  array, 
and  computing  the  area  under  this  curve  in  subroutine  AREAF. 
AREAF  is  also  used  to  find  the  area,  TIAREA,  under  the  beta 
curve  for  the  most  significant  data  Nl,  Cl  (this  area  should 
be  close  to  or  exactly  equal  to  1.0  and  is  computed  and  used 
as  a  self  adjusting  calculation  to  account  for  digital  round¬ 
off  error).  K  is  then  set  equal  to  CMAREA  divided  try-  TIAREA. 
The  output  of  this  subroutine  is  the  scalar  significance 
factor,  K. 
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PROGRAM  OPERATION 


This  Section  describes  the  inputs  necessary  to  run  the  program 
RELY  as  well  as  a  description  of  the  output  variables.  To  allow  for 
the  desireability  of  analyzing  several  different  system  and/or  flight 
data  sets  for  a  given  set  of  component  data,  the  inputs  have  been  so 
arranged  to  facilitate  these  additions  with  a  minimum  of  card  entries. 

One  may  note,  further,  that  since  the  bulk  of  the  calculations  are 
concerned  with  component  data,  additional  results  depending  on  varying 
system  or  flight  data  sets  require  only  a  small  investment  of  computer 
running  time  for  the  added  evaluation. 

To  operate  the  program  RELY  the  probability  equation  for  the 
circuit  must  be  inserted  into  the  program  between  the  COMMENT  state¬ 
ments  MATH  M)DEL- BEGIN  and  MATH  MODEL- END  (see  program  listing).  In 
doing  this  one  should  number  each  component  type  in  any  convenient 
manner  so  that  a  correspondence  exists  between  the  component  test  results 
and  the  same  component  in  the  circuit.  An  example  for  a  simple  circuit 
is  contained  in  the  sample  inputs.  In  addition  values  for  the  variables 
listed  below  must  be  supplied  in  the  format  indicated,  as  well  as  the 
actual  component,  system  and  flight  input  data. 

VARIABLE  DESCRIPTION  FORMAT 

T  Time  between  testing  dates  (months).  12 

NRS  Number  of  times  the  circuit  equation  I 6 

is  to  be  evaluated  (NRS>l)  i.e.  number 
of  Monte  Carlo  simulations,  limited  only  • 

by  format  size  and  computer  running  time. 
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VARIABLE 

DESCRIPTION 

FORMAT 

DELTAP 

(A  P)  Desired  increment  on  P- values 
in  setting  up  Beta- distributions  and 
reliability  distribution  (.01<AP  <.5>) 

F6.U 

NCOMP 

Number  of  distinct  component  types  in 
the  circuit  (l<  NCOMP  <  100 ) 

Hi 

NSYDT 

Number  of  system  data  sets  to  be  used 
with  a  given  set  of  component  data 
(l<  NSYDT  <  100) 

Hi 

NFLDT 

Number  of  flight  data  sets  to  be  used 
in  conjunction  with  the  given  component 
data  set  (l  <  NFLDT  <  100 ) 

Ili 

Values  for  these  six  variables  form  the  first  card  of  the  data 
deck.  This  card  is  followed  by:  the  component  data  set,  which  contains 
one  card  for  each  component  type,  the  laboratory  system  data  set  and 
finally  the  flight  data  set.  Each  card  of  each  set  has  the  same  format 
(hl6).  The  four  values  on  these  cards  are  N.  ,C.  ,N.  ,C.  in  that 

10*  io7  li7  n 

order.  The  first  card  of  the  component  data  set  must  contain  the  data 
corresponding  to  the  first  component  as  labelled  in  the  circuit  equation, 
the  second  card  contains  data  for  component  2,  and  so  forth  for  all 
component  data. 

All  data  cards  of  the  system  and  flight  data  sets  are  set  up 
exactly  like  those  of  the  component  data  set,  the  order  of  cards  within 
each  set  being  entirely  arbitrary.  All  combinations  of  system  and  flight 
data,  will  be  evaluated  when  NSYDT>1  and  NFLDT>1,  all  cases  of  flight 
data  being  combined  in  turn  with  each  case  of  system  data.  The  program 
may  be  run  without  system  and/or  flight  data.  In  this  case  NSYDT  and/ 
or  NFLDT  must  be  set  equal  to  1,  and  the  appropriate  data  set  must 
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consist  of  one  blank  card. 


For  a  description  of  the  output  of  program  RELY,  refer  to  the 
sample  case,  page  25. 
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SAMPLE  CASE 


As  a  sample  case  consider  the  simple  circuit  configuration  below. 

It  consists  of  seven  components  but  only  three  distinct  component  types: 
type  1,  type  2,  type  3. 


The  probability  equation  for  this  circuit  is 

R  =  (1  -  (1  -  P(l))3)  (1  -  (1  -  P( 2 ) ) 2 )  (1  -  (1  -  P(3))2) 
s 

where  P(i)  is  the  probability  of  component  type  i  functioning  properly. 

Typical  input  data  for  this  circuit  is  presented  in  the  table  on 
the  next  page.  It  will  be  noted,  for  example,  that  components  of  type  2 
worked  satisfactorily  b2  times  out  of  ljlj.  trials  during  the  first  round 
of  testing  (Tq)  but  twenty  four  months  later  (T  ),  worked  only  38  times 
in  U;  trials. 

In  order  to  compute  an  estimate  of  overall  system  reliability  for 
this  circuit  based  on  the  data  compiled  in  the  following  table  using 
program  RELY  it  is  necessary  to: 

1.  Insert  the  circuit  equation  into  the  program  between  the 
comment  cards  MATH  MODEL  -  BEGIN  and  MATH  MODEL  -  END  (See 
Listing  which  contains  this  example).  * 
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FIRST  ROUND  OF 

TESTING  (TIME  T  ) 

0 

SECOND  ROUND  OF 
TESTING  (TIME  T  ) 

i 

Number 
of  trials 

Number  of 
successes 

Number 
of  trials 

Number  of 
successes 

COMPONENTS: 

TYPE  1 

66 

Si 

66 

Sh 

TYPE  2 

hh 

h2 

hh 

38 

TYPE  3 

hh 

ho 

hh 

35 

LABORATORY 

SYSTEM  DATA: 

FIRST  CASE 

22 

22 

22 

22 

SECOND  CASE 

22 

22 

22 

18 

FIELD 

SYSTEM  DATA 
(only  one  case  was 
used  for  this 
particular  run) 

k 

1 

h 

h 

3 

Input  Data 
for 

Sample  Case 


2.  Set  up  a  data  deck  consisting  of  (refer  to  fig.  4): 


a.  One  card  containing  6  input  parameters. 


Columns  Value 


T 

1-2 

24 

NRS 

3-8 

10000 

DELTAP 

9-14 

.01 

NC0MP 

15-18 

3 

NSYDT 

19-22 

2 

NFLDT 

23-26 

1 

It  is  assumed  twenty  four 
months  elapsed  between  the 
two  rounds  of  testing.  (T 
used  only  for  printout ) . 

Number  of  simulations. 

Here  the  circuit  equation  will 
be  evaluated  10,000  times  in 
determining  Nceq,  Cceq. 

AP-  Increment  on  P 
This  will  cause  the  Beta- 
distributions  and  reliability 
histogram  to  be  formed  in 
steps  of  .01  along  the 
horizontol-axis,  P. 

Number  of  distinct  component 
types. 

There  are  three  distinct  types 
of  components  in  this  circuit 
even  though  there  is  a  total 
of  seven  components. 

Number  of  cards  in  the  system 
data  set. 

Since  there  are  two  cases  of 
laboratory  system  data,  there 
must  be  two  data  cards  in  the 
system  data  set,  as  governed 
by  this  2. 

Number  of  cards  in  flight  data 
set. 

For  this  particular  run  there 
is  only  one  case  of  flight 
data  and  consequently  there  is 
only  card  in  the  flight  data 
set.  ' 
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FIGURE  h 


DATA  DECK  FOR  SAMPLE  CASE 


b.  Component  Data  Set  -  consists  of  three  data  cards,  one 
for  each  of  the  three  distinct  component  types  in  the  sample  circuit. 

Card  1  contains  the  data  for  the  component  type  labeled#l  in  the  cir¬ 
cuit  diagram  and  circuit  equation;  card  2  contains  the  data  for  compon¬ 
ent  type  2;  and  card  3  contains  data  for  component  type  3.  This 
correspondence  between  component  numbering  in  the  circuit  diagram  and 
circuit  equation  and  the  ordering  of  cards  in  the  data  deck  is  essential, 
however  the  initial  assignment  of  numbers  to  the  different  component 
types  in  the  circuit  is  completely  arbitrary.  Each  card  in  this  data 
set  contains  the  appropriate  four  values  taken  from  the  table  of  test 
results  above.  For  example,  card  1  contains  the  number  of  trials  66 

and  corresponding  number  of  successes,  57,  for  component  1  from  the 
first  round  of  testing  at  time  Tq,  followed  by  the  test  results  from 
round  two,  time  T  :  66,  5b.  The  numbers  are  typed  in  the  first  2b 
columns  (6  columns  per  number,  right-adjusted  format-lil6). 

c.  System  Data  Set  -  consists  of  two  data  cards:  one  for  each 

of  the  two  cases  of  laboratory  system  tests.  For  example,  the  first 

? 

card  in  this  set  contains  the  results  of  the  first  case  compiled  during 

the  first  round  of  tests,  time  T  -  22  trials,  22  successes,  followed 

^  0  '  * 

by  those  of  the  second  round  of  tests,  time  T  -  22  trials,  22  successes. 
The  format  of  the  cards  in  this  set  is  identical  to  that  of  the  cards  in 
the  component  data  set. 

d.  Flight  Data  Set  -  this  data  is  set  up  exactly  like  the 
previous  two  sets  except  that  the  data  are  the  results  of  testing  the 


2k 


system  in  the  field. 


The  output  of  program  RELY  as  exhibited  in  the  output  listing  of 

the  sample  case,  Figure  5,  consists  of  a  listing  of  the  input  values 

in  tabular  form  which  shows  the  order  in  which  the  system  and  flight 

data  are  evaluated.  The  system  reliability  histogram  showing  the 

distribution  of  the  Rs  values  calculated  from  the  circuit  equation  is 

exhibited  next.  The  mean  and  variance  are  also  printed  from  which 

N  and  C  are  computed.  Finally  the  calculated  values  for  all 
ceq  cgq 

the  A.  and  intermediate  variables  as  well  P,  , ,  Mu  and  Variance 
leq  best* 

for  all  cases  of  system  and  flight  data  combinations  are  listed.  The 
printout  of  the  histogram  begins  with  the  line  preceeding  the  first 
appearance  of  a  non-zero  element  in  the  distribution.  For  example,  in 
the  output  of  the  sample  case,  the  smallest  R  value  computed  in  the 
10,000  simulations  fell  in  the  range  .77<R  <.78,  so  the  first  line 
printed  was  for  P  =  .77. 


I 
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